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Abstract. Recently many new classes of integrable systems in n dimensions occurring in 
classical and quantum mechanics have been shown to admit a functionally independent set 
of 2n — 1 symmetries polynomial in the canonical momenta, so that they are in fact su- 
perintegrable. These newly discovered systems are all separable in some coordinate system 
and, typically, they depend on one or more parameters in such a way that the system is 
superintegrable exactly when some of the parameters are rational numbers. Most of the con- 
structions to date are for n = 2 but cases where n > 2 arc multiplying rapidly. In this article 
we organize a large class of such systems, many new, and emphasize the underlying mecha- 
nisms which enable this phenomena to occur and to prove superintegrability. In addition to 
proofs of classical superintegrability we show that the 2D caged anisotropic oscillator and 
a Stackel transformed version on the 2-sheet hyperboloid are quantum superintegrable for all 
rational relative frequencies, and that a deformed 2D Kepler-Coulomb system is quantum 
superintegrable for all rational values of a parameter k in the potential. 
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1 Introduction 



We define an n-dimensional classical superintegrable system to be an integrable Hamiltonian 
system that not only possesses n mutually Poisson - commuting constants of the motion, but in 
addition, the Hamiltonian Poisson-commutes with 2n — 1 functions on the phase space that are 
globally defined and polynomial in the momenta. Similarly, we define a quantum superintegrable 
system to be a quantum Hamiltonian which is one of a set of n independent mutually commuting 
differential operators, and that commutes with a set of 2n — 1 independent differential opera- 
tors of finite order. We restrict to classical systems of the form T-l = "^^j^i g^-^PiPj + V and 
corresponding quantum systems H = A„ -|- V. These systems, including the classical Kepler and 
anisotropic oscillator problems and the quantum anisotropic oscillator and hydrogen atom have 
great historical importance, due to their remarkable properties, [H [2]. The order of a classical 
superintegrable system is the maximum order of the generating constants of the motion (with 
the Hamiltonian excluded) as a polynomial in the momenta, and the maximum order of the 
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quantum symmetries as differential operators. Systems of 2nd order have been well studied 
and there is now a structure and classification theory [21 HI O El 13 IB]- For 3rd and higher 
order superintegrable systems much less is known. In particular there have been relatively few 
examples and there is almost no structure theory. However, within the last three years there has 
been a dramatic increase in discovery of new families of possible higher order superintegrable 
classical and quantum systems [HI [101 ttH US [HI [Ml [HI [16] . The authors and collaborators have 
developed methods for verifying superintegrability of these proposed systems, [T71 [18l [191 120] . 
In the cited papers the emphasis was on particular systems of special importance and recent 
interest. Here, however, the emphasis is on the methods themselves. 

In Section [2] we review a method for constructing classical constants of the motion of all 
orders for n-dimensional Hamiltonians that admit a separation of variables in some orthogonal 
coordinate system. Then we apply it to the case n = 2 to show how to derive superintegrable 
systems for all values of a rational parameter k in the potential. Many of our examples are new. 
In Section [3] we review our method for establishing a canonical form for quantum symmetry oper- 
ators of all orders for 2-dimensional Schrodinger operators such that the Schrodinger eigenvalue 
equation admits a separation of variables in some orthogonal coordinate system. Then we apply 
this method to establish the quantum superintegrability of the caged anisotropic oscillator for all 
frequencies that are rationally related, and for a Stackel transformed version of this system on 
the 2-sheet hyperboloid. We give a second proof for the caged oscillator in all dimensions n that 
relies on recurrence relations for associated Laguerre polynomials. We also apply the canonical 
equations to establish the quantum superintegrability of a deformed Kepler-Coulomb system 
for all rational values of a parameter k in the potential. 

1.1 The construction tool for classical systems 

There are far more verified superintegrable Hamiltonian systems in classical mechanics than was 
the case 3 years ago. The principal method for constructing and verifying these new systems 
requires that the system is already integrable, in particular, that it admits a separation of 
variables in some coordinate system. For a Hamiltonian system in 2?i-dimensional phase space 
the separation gives us n second order constants of the motion in involution. In this paper we 
first review a general procedure, essentially the construction of action angle variables, which 
yields an additional n — 1 constants, such that the set of 2n — 1 is functionally independent. 
This is of little interest unless it is possible to extract n new constants of the motion that are 
polynomial in the momenta, so that the system is superintegrable. We will show how this can be 
done in many cases. We start first with the construction of action angle variables for Hamiltonian 
systems in n dimensions, and later specialize to the case n = 2 to verify superintegrability. 

Consider a classical system in n variables on a complex Riemannian manifold that admits 
separation of variables in orthogonal separable coordinates xi, . . . Then there is an n x n 
Stackel matrix 



S = {Sij{xi)) 



such that $ = det S ^ and the Hamiltonian is 



n 



n 



n = Ci = '^Tu [pf + Vi{xi 




i=l 



i=l 



where V = Y17=i Vi{xi), T is the matrix inverse to S: 



n 



n 






i=i 
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and 6ik is the Kronecker delta. Here, we must require H^^^Tij / 0. We define the quadratic 
constants of the motion C^, k = 1, . . . ,n hy 

n 

Ck ='Y2'^ki{pl + Vi), k = l,..., n, 

i=l 

or 

n 

Pi + = X] "^^i^i' 1 < i < n. (2) 

As is well known, 

{Cj,Ck} = 0, l<j,k<n. 

Here, 

n 

{^(x, p), i3(x, p)} = ^{d,A dp^B - dp^A diB). 

i=l 

Furthermore, by differentiating identity ([1]) with respect to we obtain 

n 

dhTit = - ^ TihS'f^jTji, I < h,i,i < n, 
where S'f^j = dhShj. 

Now we define nonzero functions Mi:j{xj,pj,Ci, . . . ,Cn) on the manifold by the requirement 
{Mkj,Ce} = TijSjk, I <k,j,e < 71. 
It is straightforward to check that these conditions are equivalent to the differential equations 



n. 



We can use the equalities ([2]) to consider M^j either as a function of xj alone, so that -^M^j = 
Sjk/2pj or as a function of pj alone. In this paper we will take the former point of view. 
Now define functions 

n 

Cg = Y,Mgj, l<q<n. 
Then we have 

n 

{CqjCl} = ^ TijSjq = 5lq. (3) 

i=i 

This shows that the 2n — 1 functions 
Ti. = £i, £21 • • • 1 C,n, -^21 • • • 5 C,n, 
are constants of the motion and, due to relations ([3]), they are functionally independent. 
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Now let's consider how this construction works in n = 2 dimensions. By replacing each 
separable coordinate by a suitable function of itself and the constants of the motion by suitable 
linear combinations of themselves, if necessary, e.g. [21], we can always assume that the Stackel 
matrix and its inverse are of the form 



S 



fi 

f2 



1 



-1 



T 



1 



/1 + /2 



1 1 

/2 -/l 



where fj is a function of the variable Xj alone. The constants of the motion Ci = % and C2 
are given to us via variable separation. We want to compute a new constant of the motion C2 
functionally independent of £1, £2- Setting M21 = M, M22 = —N, we see that 



2pi-^M = 1, 
ax I 



2p2-^N 

dX2 



1, 



(4) 



from which we can determine M, N. Then £2 = — is the constant of the motion that we 
seek. 

The treatment of subgroup separable superintegrable systems in n dimensions, |18) . is also 
a special case of the above construction. Suppose the Stackel matrix takes the form 





( 1 


-/i 


, ••• 




\ 









1 


-/2 









s = 

























••• 


1 


-fn-1 






^ 





••• 





1 / 




where fi = 


fi{Xi) 


, i = 


1, . . . , n. The inverse matrix in this case is 




( 1 


fi 


/1/2 ••• 




/1/2 • • • fn-1 \ 









1 


/2 72/3 




/2 • • • /n 


-1 




T = 





























1 


fn-1 








V 











1 


/ 





which leads exactly to the construction found in 



1.2 Application of the construction for n = 2 

As we have seen, for n = 2 and separable coordinates ui = x, U2 = y we have 

1 



£2 



fi{x) + f2{y) 
f2{y) / 2 



[Px+Py + Vl{x) +V2{y)), 



fi{x) + f2{y) 



p^ + vi{x)) 



fi{x) 



fi{x) + f2{y) 



{pl + v2{y)) ■ 



We will present strategies for determining functions /i, /2, vi, vi such that there exists a 3rd 
constant of the motion, polynomial in the momenta. The constant of the motion C2 = M — N 
constructed by solving equations (jlj) is usually not a polynomial in the momenta, hence not 
directly useful in verifying superintegrability. We describe a procedure for obtaining a polynomial 
constant from M — N, based on the observation that the integrals 



M 



1 



dxi 



2 J Vfi^ + C2-V1 



N 



dxo 



V/2^ - £2 



V2 
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can often be expressed in terms of multiples of the inverse hyperbolic sine or cosine (or the 
ordinary inverse sine or cosine), and the hyperbolic sine and cosine satisfy addition formulas. 
Thus we will search for functions fj, vj such that M and possess this property. There 
is a larger class of prototypes for this construction, namely the second order superintegrable 
systems. These have already been classified for 2-dimensional constant curvature spaces [1] and, 
due to the fact that every superintegrable system on a Riemannian or pseudo-Riemannian space 
in two dimensions is Stackel equivalent to a constant curvature superintegrable system |22| 123]. 
this list includes all cases. We will typically start our construction with one of these second 
order systems and add parameters to get a family of higher order superintegrable systems. 

A basic observation is that to get inverse trig functions for the integrals M, we can choose 
the potential functions f{z), v{z) from the list 

.9 „ ^ , N A + B sin pz 

Zi{z) = Az^ + Bz + C, Z2{z) = 5 —, 

cos^ py 

Zsiz) = — ^ + — f— , Z^{z) = Ae^'P' + Be'P', 
cos^ pz sm pz 

„ / ^ .9 B „ / s A B 
Z5{z) = Az^ + -^, Zg{z) = - + -^, 

sometimes restricting the parameters to special cases. Many of these cases actually occur in 
the lists [5 so we are guaranteed that it will be possible to construct at least one second order 
polynomial constant of the motion from such a selection. However, some of the cases lead only 
to higher order constants. 

Cartesian type systems 

For simplicity, we begin with Cartesian coordinates in flat space. The systems of this type are 
related to oscillators and are associated with functions Zj for j = 1, 5, 6. We give two examples. 

[El]. For our first construction we give yet another verification that the extended caged har- 
monic oscillator is classically superintegrable. We modify the second order superintegrable 
system [El], [4], by looking at the potential 

rf-Z yZ 

with £2 = + wfx^ + (For uji = 0J2 this is just [El].) It corresponds to the choice of 
a function of the form Z5 for each of vi, V2- Evaluating the integrals, we obtain the solutions 

M{x,p^) = —A, smhA = ^7=1= ^ 



4a;i y/C'l - Aco'iP 

N[y,Py) = -j — smhB - 



4W2 ' V(?^- £2)2-4^^^7' 

to within arbitrary additive constants. We choose these constants so that M, N are proportional 
to inverse hyperbolic sines. Then, due to the formula cosh^ u — sinh^ u = 1, we can use the 
identities ([2]) to compute cosh^ and coshi3: 

2ujixp.^ _ ^ 2uj2ypy 



cosh^= — , ^ =, coshS 



hat iOi/ijJ2 = 
Then coi = pco, UJ2 = qoj and 



Now suppose that U}i/uj2 = k is rational, i.e. k = ^ where p, q are relatively prime integers. 



sinh(-4ip^a;[M - A^]) = sinh(^^ + pB), cosh.{-Aipquj[M - A^]) = cosh(g^ + pB), 
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are both constants of the motion. Each of these will lead to a polynomial constant of the motion. 
Indeed, we can use the relations 

(coshx lb sinhx)" = coshnx it sinhnx, 
cosh(2; + y) = cosh x cosh y + sinh x sinh y, 
sinh(x + y) = cosh x sinh y + sinh x cosh y, 

coshnx = ( I sinh^-' x cosh""^-' x, 

[(n=l)/2] 



sinh ?ix = sinh X I „. - 1 sinh^-'"^ x cosh"~^-'~^ x. 



recursively to express each constant as a polynomial in cosh^, sinh^, cosh^S, sinh;B. Then, 
writing each constant as a single fraction with denominator of the form 

we see that the numerator is a polynomial constant of the motion. Note that, by construction, 
both sv'Qh.{—Aipqu[M — N]) and cos\i{—Aipqu][M —N]) will have nonzero Poisson brackets with C21 
hence they are each functionally independent of 7^, £2- Since each of our polynomial constants 
of the motion differs from these by a factor that is a function of £2 alone, each polynomial 
constant of the motion is also functionally independent oil-L, £2- Similar remarks apply to all 
of our examples. 

[E2] . There are several proofs of superintegrability for this system, but we add another. Here, we 
make the choice vi = Zi, V2 = Z^, corresponding to the second order superintegrable system [E2] 
in [1]. The potential is 

TA 2 2,2 2, , 

v=uj^x + ujoy + Qx H — ^, 

where £2 = pt + uifx'^ + ax and the system is second order superintegrable for the case coi = uji. 
Applying our method we obtain 

M{x,p^) = —A, smhA = ^ - , cosh^- ^ 



Niy,Py) = — 'B, sinhi3 = , coshi3 



4w2 ^{u- C2f -Auj'ip y/{n - C2? - 4a;i/3 

Thus if U1/2LJ2 is rational we obtain a constant of the motion which is polynomial in the 
momenta. 

Polar type systems 

Next we look at flat space systems that separate in polar coordinates. The Hamiltonian is of 
the form 

n=pl + ^ {pI + f{r) + g{e)) = {pj, + pj + vi (R) + V2{e)) , 
x = i? = lnr, y = 9, /i = e^^, /2 = 0. 
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Cases for which the whole process works can now be evaluated. Possible choices of / and g are 



a 



{l)f{r)=ar', g{e) = Z^{k9), (2) /(r) = -, g{e) = Z-sikO) 



r 
a 



{2,)f{r) = ar\ g{9) = Z^{k9), (4) /(r) = -, g{e) = Z^{ke), 

^5)f^r) = ar^ g{9) = Z2{ke), (6) /(r) = ^, g{e) = Z^ikO). 

In each case if p is rational there is an extra constant of the motion that is polynomial in the 
canonical momenta. 

[El]. Case (1) is system [El] for k = \. For general k this is the TTW system [12], which we 
have shown to be supperintegrable for k rational. 

Case (2). This case is not quadratic superintegrable, but as shown in [19] it is superintegrable 
for all rational k. 

[E8]. For our next example we take Case (3) where z = x + iy: 

for arbitrary k. (li k = 2 this is the nondegenerate superintegrable system [E8] listed in [5].) In 
polar coordinates, with variables r = and z = e^^*^. Then we have 

-£2 =pU ^ie'''' + le^"', n = e"2«(p| - £2 + 4ae^«) . 
Our method yields 

Nie,pe) = -^B, sinh^ = + cosh^ = J^Ei^^,-^^o 

M(i?,p«) = -^A, sinh^ = i-^^^^"''' - H) ^^^^^ ^ 2^pK^^^_2« 

This system is superintegrable for all rational k. 
[E17]. Taking Case (4) we have, for z = x + iy: 



J zz z^~^^ 2:^/2+1 
(For k = 1 this is the superintegrable system [E17] in [3].) Then we have 

H = e-2« +pI + + fie-^'^' + 76^*^^) , (5) 
-£2 =pI + fie-^'""' + 7e-^'=^ n = e-'^{p%, - £2 + ae«) . 
Applying our procedure we find the functions 

Nie,pe) = ^B, sinh^ = (l^^lM, cosh^ = _^E^e-^ 
^ ' k^2 ' V4/3£2-72 

M{R,pn) = ^A, smhA=t^£^L^±^, cosh^ = ^2^£=e-«. 
This demonstrates superintegrability for all rational k. 

[E16]. Case (5) corresponds to [E16] for k = 1, and our method shows that it is superintegrable 
for all rational k. 

Case (6). This case is not quadratic superintegrable, but it is superintegrable for all rational k. 
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Spherical type systems 



These are systems that separate in spherical type coordinates on the complex 2-sphere. The 
Hamiltonian is of the form 



x = y = iIj, /i = 0, /2 = 



1 



cosh 



Embedded in complex Euclidean 3-space with Cartesian coordinates, such 2-sphere systems can 
be written in the form 



where 



Js = SiPs2 - S2Psi , Jl = S2Ps 



S3Ps 



J2 = SsPsi - SlPs3, 



sf + + s'i 



[S9]. Here, we have the case vi{ip) = and V2{ip) is a special case of Z3. 



V 



a B -y 



It is convenient to choose spherical coordinates 

cos tp 



Si 



S2 = simp cosh^, S3 = tanh^. 



sinh ip ' 

In terms of these coordinates the Hamiltonian has the form 



Ti = cosh ^p 



a 



+ 



/3 



7 



cos^ ip sin^ ip sinh^ ip 



a 



+ 



/3 



(n ' 9 1.95 

^ cos^ ip sm (f 



H = cosh iIj 



PI+C2+ . \ 
smh yj 



system [S9] in [^. We extend this Hamiltonian via the replacement ip ^ kip and proceed with 
our method. Thus 



T-L = cosh ^ 



a 



a 



+ 



/3 



7 



cos^ kip ' sin^ k(p sinh^ V 

/5 



with Ti expressed as above. The functions that determine the extra constant are 



M{^,p^) 
sinh^ = 

N{iJ,p^) 

sinhi3 = 



--A, 



Aky/C2 

i{C2 cos{2kip) — a + /3) 
V^(£2 - 4q/3' 



cosh^ 



sm{2k(p)pi^ 



i{C2 cosh(2V^) + 7 - -?^) 

y/{n + C2-jy+4.C27 



coshB 



^{£2 - a - /3)2 - 4a/3 
i sinh(2^/')p^ 

" y/{n+C2 --fy+4C2j' 



Thus this system is superintegrable for all rational k. 
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[S7]. This system corresponds to vi{(p) = Z2 and V2{ip) a variant of Z3. The system is second 
order superintegrable: 

V4 + 4 sWs'i + s'i si' 

Choosing the coordinates ip and ip we find 

u u2 , 2 , 2 , sinhy; cosy? 7 

7^ = cosh iplp^+p +a TTT+P— ^ ^ . 2 ' 

V cosh ip sm (f sm ip J 

^2=P^ + P— ^ h — o— , = cosh V + ^2 + a -^j— . 

sm ip sm 99 \ cosh ip J 

We make the transformation ip — t- and obtain 

^•2,(1 2 sinh „ cos kip 1 \ 

■H = cosh ip\v^^V^^ OL r^+/3— h . 2 , ' 

V cosh ^ sm kip sm fcy? / 

2/1 cos kip 7 

^2 = + ^ . 2r + ■ 2. ' 

sm kip sm /ec^ 

with as before. The functions that determine the extra constants are 
\/ L2k 

i{£2Cos{kip) + 13) ^^^^ ^_ 2^/C^sm{kip)p^ 



sinh A = — : , cosh A — , „ , 

V L2 

. , „ 2L2 sinh ■0 + a , „ 2i cosh V'p^ 

smh Z3 = — — , cosh B — 



V-a2 + 4£2 _ 4^2^ ' ^-a2 + - 4^2^ ' 

Thus this system is superintegrable for all rational k. 

[S4]. Here vxiip) = Z4 and V2{ip) is a variant of Z2. This is another system on the sphere that 
is second order superintegrable and separates in polar coordinates: 



{si-is2)^ y^s'i + s'i {si-is2)^/sl + s{ 
In terms of angular coordinates ip, ip the Hamiltonian is 

n = cosh^ ip (4 +pI + ae^^^^ + + ^^^) • 
After the substitution ip ^ kip we have 

C2=pI + ae'^'^ + y'^^ U = cosh^ ip [pi + C2 + ■ 
The functions that determine the extra constants are 

Miip^p,) = -^A, smh^ = cosh^ = 

iV(0,p,) = ^^, sinhg= £2 sinh 0-/3 ^ ^^^^^ _ 2.cosh0p, 



' ^4£i + 4£2-/32' V4£i + 4£2-/32' 

so this systems is also superintegrable for all rational k. 
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[S2]. Here vi{(p) = Z4 and V2{tp) is a variant of Z3: 
_a /3 7(si + is2) 



(si-is2)2 {si-is2y 
which is is second order superintegrable. After the substitution ip — )• we obtain the system 

n = COsh^ (4 ^ ^2 ^ _|_ ^ f^^2^k^ ^ ^g4^/c^^ ^ (g) 

C2=pI + Pe^'"'^ + 76^'^=^ , n = cosh^ ^j.^^ — ^) • 

Applying our procedure we find 

sinh^ = ^(^2Cosh(2V.)-a + 7^) ^^^^^ _ V^sinh(2V.)p^ 



^(£2 - a - 7^)2 - 4a?^ ^(£2 -a-nf- AaH 

Thus the system is superintegrable for all rational k. 

1.3 Horospherical systems 

In terms of horospherical coordinates on the complex sphere we can construct the Hamiltonian 

n = y'^ {pI +pI + + + a + px). 

If uji/uj2 is rational then this system is superintegrable. However, there is no need to go into 
much detail for the construction because a Stackel transform, essentially multiplication by 1/y^ , 
takes this system to the flat space system generalizing [E2] and with the same symmetry algebra. 
There is a second Hamiltonian which separates on the complex 2 sphere 

TT 2 ( 2 , 2 , , Q , 2 2, 2 2\ 

H = y \p^+Py + ^ + 13 + ujix +uj2y )■ 

It is superintegrable for uji/ijJ2 rational, as follows from the Stackel transform l/y^ from the 
sphere to flat space. If oji = 0J2 then we obtain the system [S2]. We note that each of the 
potentials associated with horospherical coordinates can quite generally be written in terms 
of si, S2j 53 coordinates using the relations 

~^ 2,2,2 1 

y= — , X = — , Si + S2 + S3 = l. 

Si - IS2 Si - IS2 

1.4 Generic systems 

These are systems of the form 

where j, k can independently take the values 2, 3, 4 and Zj^ Zj depend on distinct parameters. 
Superintegrability is possible because of the integrals 



/ 



dx 1 /2Ce-2^^'^' + S 

arcsin 
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dx 1 / -2Csmkx + B \ 

— — arcsin = 
kVC W4.C^ + 4:CA + B^J 



A+B sin kx 



cos 



kx 



I 



dx 1 I C cos 2kx + A-B 

: arcsin 



_A^B^^ 2kVC \y/{A + B + Cy-4AB 

cos^ kx sin^ kx ^ ' 



Consider an example of this last type of system: 

-X/ — \ /'„2 _ „2 I A 2ikx , tj 4ikx _ „p2iqy _ hp'^iqy] 

The equation % = E admits a separation constant 

pi + {A- EAy^'^"" + {B- EB')e^'''^ = + {a - Ea')e^''^y + {b - Eh')e^''^y = C. 
As a consequence of this we can find functions M{x,px) and N{y,py) where 

N 1 ■ f 2Ae-'^'^'= + b\ 

M{x,Px) = ^ arcsm — — , 

AkVC [ y/^LA - B^ J 

where A = A- EA', B = B- EB' and 

where a = a — Ea' , b = b — Eb' . We see that if | is rational then we can generate an extra 
constant which is polynomial in the momenta. 

The superintegrable systems that have involved the rational functions Zi, Z5 and Zq work 
because of the integrals 

dx 1 . , f 2Ax + C 

arcsmh 



and 



/ 



y/Ax^ + Bx + C y/A 



dx 1 / 2B + Ax \ 

, = ■= arcsmh — = 

+ Bx + C Vb \xVA^^^4BCJ 



2 Quantum superintegrability 

2.1 The canonical form for a symmetry operator 

We give a brief review of the construction of the canonical form for a symmetry operator |19j . 
Consider a Schrodinger equation on a 2D real or complex Riemannian manifold with Laplace- 
Beltrami operator A2 and potential V: 

= (A2 + V)^ = E^ 

that also admits an orthogonal separation of variables. If {ni,U2} is the orthogonal separable 
coordinate system the corresponding Schrodinger operator has the form 

H = L, = A2 + V{m,u,) = + 9l + vi[u,) + V2{U2)) . (7) 
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and, due to the separability, there is the second-order symmetry operator 
i.e., [L2,H] = 0, and the operator identities 

h{ui)H + L2 = + Vi{ui), f2{u2)H -L2 = + V2{U2). 

We look for a partial differential operator L{H, L2,ui,U2) that satisfies 
[H,L] = 0. 

We require that the symmetry operator take the standard form 

L = Y,{^''Hui,U2)d^,u2+B^^Hui,u2)du,+C^'''{ui,U2)du,+D^^''{ui,U2))W4. (8) 

We have shown that we can write 

L{H,L2,ui,U2) = A{ui,U2,H, 12)812 + B{ui,U2,H,L2)di 

+ na, H, 12)82 + D{ui,U2,H, L2), (9) 

and consider L as an at most second-order order differential operator in ui, U2 that is analytic 

in the parameters H, L2. Then the conditions for a symmetry can be written in the compact 
form 

-^uiui ~\~ -Au2U2 2-6^2 ~l" 2Cti]^ = 0, (10) 

Bu,u, + - '2Au2V2 + 2Du, - Av'2 + {2A^j2 + Af^)H - 2A^,L2 = 0, (11) 

C„i«i + Cu2U2 - 2Au,vi + 2Du2 - Av[ + (2A, Ji + Af[)H + 2Au,L2 = 0, (12) 
Du-^u-, + Du2U2 - 2Bu^vi - 2Cu2V2 - Bv[ - Cv'^ 

+ {2BuJi + 2Cuj2 + Bf[ + Cf^)H + {2Bu, - 2Cu2)L2 = 0. (13) 

We can further simplify this system by noting that there are two functions F{ui,U2, H, L2), 
G{ui,U2, H, L2) such that pO|) is satisfied by 

A = F, B = -^d2F-diG, C = -^diF + d2G. (14) 

Then the integrability condition for (llip . ()12p is (with the shorthand djF = Fj, dj^F = Fj£, 
etc., for F and G), 

2G1222 + ^i^2222 + 2F22(^;2 " + L2) + 8^2(1;^ - f^H) + F{v'^ - f^H) 

- 2G1112 + ^Fnn + 2Fn{vi - fiH - L2) + 2,Fi{v[ - f[H) + F{v'l - f{H), (15) 
and equation (fT3|) becomes 

^^1112 + 2Fi2(?;i - fiH) + - f'^H) + Igiiii + 2Gii(?;i - /lif - L2) 

+ Gi{v[ - f[H) = -^Fl222 - 2Fi2(7;2 - f2H) - F2{v[ - f[H) + 1^2222 

+ 2G22(^;2 - f2H + L2) + G2{v'2 - f^H). (16) 
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Here, any solution of (jlSp . (jl6p with A, B, C not identically corresponds to a symmetry 
operator that does not commute with L2, hence is algebraically independent of the symmet- 
ries H, L2. (Informally, this follows from the construction and uniqueness of the canonical form 
of a symmetry operator. The operators L = g{H, L2) algebraically dependent on H and L2 are 
exactly those such that A = B = C = Q, D = g{H, L2). A formal proof is technical.) Note also 
that solutions of the canonical equations, with H, L2 treated as parameters, must be interpreted 
in the form ([8|) with H and L2 on the right, to get the explicit symmetry operators. 

2.2 The caged anisotropic oscillator 

In [19] we used the canonical form for symmetry operators to demonstrate the quantum superin- 
tegrability of the TTW system [121ll3j for all rational k, as well as a system on the 2-hyperboloid 
of two sheets [19] . Here we give another illustration of this construction by applying it to the 2D 
caged oscillator [11]. For p = q This is the second order sup erintegr able system [El] on complex 
Euclidean space, as listed in [4]. Here, 

= {811+822 + V{ui,U2))^, (17) 

where 

V{Ui, U2) =UJ [p + q U2) + ^ + -2, 

Uf U2 

in Cartesian coordinates. We take p, q to be relatively prime positive integers. Thus 

f 1 f n 222"! 222, "2 

71 = 1, 72 = 0, vi=ujpui + ^, V2=ujqu2 + ^. 

The 2nd order symmetry operator is 

L2 = - {di + V2{u2)) , 

and we have the operator identities 

8l = -{vi{ui) + H + L2), 8l = -V2{u2) + L2. 

Based on the results of [18] for the classical case, we postulate expansions of F, G in finite 
series 

F = '^Aa,bEa,b{ui,U2), G = Ba,bEa,b{ui, U2) , Ea,b{ui, U2) = U1U2. (18) 

a,b a,b 

The sum is taken over terms of the form a = uq + m, b = bo + n, and c = 0, 1, where m, n 
are integers. The point (og, 60) could in principle be any point in M?, 

Taking coefficients with respect to the basis (I18p in each of equation (llSp and ()16p gives 
recurrence relations for these coefficients. The shifts in the indices of A and B are integers and 
so we can view this as an equation on a two-dimensional lattice with integer spacings. While 
the shifts in the indices are of integer size, we haven't required that the indices themselves be 
integers, although they will turn out to be so in our solution. The 2 recurrence relations are 
of a similar complexity, but rather than write them out separately, we will combine them into 
a matrix recurrence relation by defining 
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Figure 1. The template. Points contributing to the recurrence relation are marked with large dots (•, o). 
The large dot on the bottom center corresponds to the position (a — 1,6+ 1), the large dot at the top 
corresponds to the position (a + 4, b) and o corresponds to position (0, 0). 

We write the 2 recurrence relations in matrix form as 

= Ma^bCa^b + Ma-l^b+lCa~l,b+l + -^a,b+2C'a,fe+2 + Afa,fe+4C'a,6+4 

+ Afa+l,fe-lC'a+l,6-l + ^a+l.fc+lC'a+l.b+l + ^•^a+l,b+3C'a+l,fe+3 + ^Ia+2,bC a+2,b 
+ ^^a+3,fe+lC'a+3,6+l + ^^a+A,bC a+A,b^ 

where each Mjj- is a 2 x 2 matrix given below. 

_ / 2a;2([6+l]g- [a+l]p)([6+l]g+[a+l]p) 

\ -2u:'^{bq-ap){hq + ap) 

M -( ^ 



Ma. 



b+2 



2L2(6 + 2)(6+l) 

-2L2(6 + 2)(6 + l) 



_/ i(6 + 3)(6 + l)(62 + 66 + 4a2 + 8) 

-i(6 + 4)(6 + 2)(62+46 + 4a2 + 3) 

Ma^i,b-i=[ 2u^q%{a + l) 



M. 



a+l,b+l 




-2H{a + l){h + l) 



2(a + l)(6 + 3)(fe + 2)(6 + l) 



^"+^'''+^ ' i(a + l)(6 + 2)(62 + 46 + 4Q2 + 3) 



M, 



a+2,b 



2{H + L2){a + 2){a + l) 

-2{H + L2){a + 2){a + l) 



M -( 2(a + 3)(a + 2)(a+l)(6+l) 

M,+3,6+i - 1 (a + 2){h + l)(a2 + 4a + 4ai + 3) 

_ / -i(a + 3)(a + l)(a2 + 6o + 4ai + 8) 
Ma+4,fe - Q 1 (a + 4) (a + 2) (a2 + 4a + 4qi + 3) 

It is useful to visualize the the set of points in the lattice which enter into this recurrence 
for a given choice of (a, 6). These are represented in Fig. [TJ Although the recurrence relates 
10 distinct points, some major simplifications are immediately apparent. We say that the 
lattice point (a, h) has even parity if a + 6 is an even integer and odd parity if a + 6 is odd. Each 
recurrence relates only lattice points of the same parity. Because of this we can assume that the 
nonzero terms b will occur for points of one parity, while only the zero vector will occur for 
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points with the opposite parity. A second simplification results from the fact that the recurrence 
matrices Ma+mfi+n are of two distinct types. Either m, n are both even, in which case Ma-\.m,b+n 
is diagonal, or m, n are both odd, in which case the diagonal elements of Ma-\-m,b+n are zero. 
Another simplification follows from the observation that it is only the ratio p/q = r that matters 
in our construction. We want to demonstrate that the caged anisotropic oscillator is operator 
superintegrable for any rational r. The construction of any symmetry operator independent 
of H and L2 will suffice. By writing u = to' /2, p' = 2p, q' = 2q in H, we see that, without loss 
of generality, we can always assume that p, q are both even positive integers with a single 2 as 
their only common factor. 

If for a particular choice of p, q we can find a solution of the recurrence relations with only 
a finite number of nonzero vectors Ca^b then there will be a minimal lattice rectangle with 
vertical sides and horizontal top and bottom that encloses the corresponding lattice points. 
The top row oq of the rectangle will be the highest row in which nonzero vectors Cag^b occur. 
The bottom row oi will be the lowest row in which nonzero vectors Ca^^b occur. Similarly, 
the minimal rectangle will have right column 69 and left column 61. Now slide the template 
horizontally across the top row such that only the lowest point on the template lies in the top 
row. The recurrence gives (ao + l)Mao,fe ~ ^ ^'^^ columns b. Based on examples for the 
classical system, we expect to find solutions for the quantum system such that uq > and 
b > 0. Thus we require ^ao,6 = along the top of the minimal lattice rectangle, so that all 
vectors in the top row take the form 



a 



ao,b 





Bao,b 



Now we move the template such that the recurrence 6, second row from the bottom and on 
the left, lies on top of the lattice point (oq, bo). This leads to the requirement {boq — aop){boq + 
o-QP)Bao,bo = 0- Again, based on hints from specific examples, we postulate -Bao,6o / 0) f^o = 9) 
60 = P- Since we can assume that p, q are even, this means that all odd lattice points correspond 
to zero vectors. Next we slide the template vertically down column 5o such that only the left 
hand point on the template lies in column 6o- The recurrence gives (60 — l)o^a,6o ~ all 
rows a. Since 60 is even, we postulate that 



Ca. 



bo 





Ba,bo 



for all lattice points in column 69. 

Note that the recurrence relations preserve the following structure which we will require: 

1. There is only the zero vector at any lattice point (a, b) with a + b odd. 





2. If the row and column are both even then Ca^b 

3. If the row and column are both odd then Cab - 



Ba,b 
Aa,b 





This does not mean that all solutions have this form, only that we are searching for at least one 
such solution. 

To finish determining the size of the minimum lattice rectangle we slide the template vertically 
downward such that only the right hand point on the template lies in column bi . The recurrence 
gives (61 — l)(6i — 3)^(a, 61) = 0, 61(61 — 2)Ba^bi = for all rows a. There are several possible 
solutions that are in accordance with our assumptions, the most conservative of which is 61 = 0. 
In the following we assume only that -Bao.feo 7^ 0, the structure laid out above, and the necessary 
implications of these that follow from a step by step application of the recurrences. If we find 
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a vector at a lattice point that is not determined by the recurrences we shall assume it to be 
zero. Our aim is to find one nonzero solution with support in the minimum rectangle, not to 
classify the multiplicity of all such solutions. 

Now we carry out an iterative procedure that calculates the values of Ca^b at points in the 
lattice using only other points where the values of Cij are already known. We position the 
template such that the recurrence Ma^b, second row from the bottom and on the left (of the 
template), lies above the lattice point (oq, b), gq = q, and slide it from right to left along the top 
row. The case b = bo = p has already been considered. For the remaining cases we have 

2a;Vao(6 + l)^ao-i,b+i = '^^^{bq - aop){bq + aop)Ba,,,b + 2L2(6 + 2)(6 + l)Ba,,b+2 

+ i(6 + 4)(6 + 2) (62 + 46 + + 3)Ba,,b+A. (19) 

Only even values of b need be considered. Now we lower the template one row and again slide it 
from right to left along the row. The contribution of the lowest point on the template is and 
we find 

2u\q[b + 2] - p[ao]){q[b + 2] + p[ao])^ao-i,b+i = -2ao(6 + 4)(6 + 3)(6 + 2)Bao,b+i 

- 2L2(6 + 3)(6 + 2)Aa,-i,b+3 - ^(6 + 4)(6 + 2){b^ + 8b + lias + 8)^^0-1,6+5, (20) 

where we have replaced 6 by 6 + 1 so that again only even values of b need be considered. 
For the first step we take b = p — 2. Then equation ()19p becomes 

2a;VgA-i,P-i = -8uj'^q^iP - l)-Bq,p-2 + 2L2p{p - l)Sg,p, 

and (j20p is vacuous in this case, leaving j4g_i^p_i undetermined. Note that several of the terms 
lie outside the minimal rectangle. From these two equations we can solve for Bq p^2 in terms of 
our given i?g,p, Aq^i^p. Now we march across both rows from right to left. At each stage our two 
equations now allow us to solve uniquely for Aq_i^b+i and Bq b in terms of A's and 5's to the 
right (which have already been computed). We continue this until we reach 6 = and then stop. 
At this point the top two rows of the minimal rectangle have been determined by our choice 
of Bq^p and Aq^i^p^i. We repeat this construction for the third and fourth rows down, then the 
next two rows down, etc. The recurrence relations grow more complicated as the higher rows of 
the template give nonzero contributions. However, at each step we have two linear relations 

/ 2wVa(^ + l) -2u'^{bq-ap){bq + ap)\ ( Aa-i,b+i\ ^ 

\2u^{q[b + 2]-pa){q[b + 2]+pa) ^ V Ba,b ) '"' 

where the right hand side is expressed in terms of ^'s and -B's either above or on the same line 
but to the right of Aa-i^+i, Ba^b, hence already determined. Since the determinant of the 2x2 
matrix is nonzero over the minimal rectangle, except for the upper right corner and lower left 
corner, at all but those those two points we can compute Aa-i^+i, Bafi uniquely in terms of 
quantities already determined. This process stops with rows a = 1,2. 

Row a = 0, the bottom row, needs special attention. We position the template such that 
the recurrence M^^ lies above the lattice point (0,6), and slide it from right to left along the 
bottom row. The bottom point in the template now contributes so at each step we obtain an 
expression for B^ b in terms of quantities A, B from rows either above row or to the right of 
column 6 in row 0, hence already determined. Thus we can determine the entire bottom row, 
with the exception of the value at (0,0), the lower left hand corner. When the point M^^b on 
the template is above (0, 0) the coefficient of -Bq.o vanishes, so the value of i?o,o is irrelevant 
and we have a true condition on the remaining points under the template. However, this is 
a linear homogeneous equation in the parameters -B<j,p, Aq_i,p_i: xBq,p + r/^g_i,p_i = for 
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constants x-, V- Hence if, for example, x / we can require Bg^p = — (r//x)^q-i,p-i and satisfy 
this condition while still keeping a nonzero solution. Thus after satisfying this linear condition 
we still have at least a one parameter family of solutions, along with the arbitrary i?o,o (which 
is clearly irrelevant since it just adds a constant to the function G). However, we need to check 
those recurrences where the point Ma^b on the template slides along rows a = —1, —2, —3, —4 to 
verify that these relations are satisfied. We have already utilized the case a = —4, and for cases 
a = —3, —2, —1 it is easy to check that the relations are vacuous. 

The last issue is the left hand boundary. We have determined all vectors in the minimal 
rectangle, but we must verify that those recurrences are satisfied where the point Ma^b on the 
template slides along columns b = —1, —2, —3, —4. However again the recurrences are vacuous. 

We conclude that there is a one-parameter (at least) family of solutions to the caged aniso- 
tropic oscillator recurrence with support in the minimal rectangle. By choosing the arbitrary 
parameters to be polynomials in H, L2 we get a finite order constant of the motion. Thus the 
quantum caged anisotropic oscillator is superintegrable. We note that the canonical operator 
construction permits easy generation of explicit expressions for the defining operators in a large 
number of examples. Once the basic rectangle of nonzero solutions is determined it is easy to 
compute dozens of explicit examples via Maple and simple Gaussian elimination. The generation 
of explicit examples is easy; the proof that the method works for all orders is more challenging. 

Example. Taking p = 6 and q = 4:, we find (via Gaussian elimination) a solution of the 
recurrence with nonzero coefficients 

^1,1, ^1,3; ^1,5, ^3,1, ^3,3, ^3,5, 

Bo,2, -Bo,4, B2fi, i?2,2, -82,4, i?4,0) -64^2, -64^4, Bq^, i?6,2, Bq^^. 

We can take ^1^5 = ai and ^3^3 = 02 and then all other coefficients will depend linearly on ai 
and 02. Calculating the A, B, C and D from equations (jlip . ()12p . (jl4p we find the coefficient 
of has a factor of 2a2 + 9ai and so we set ai = co^ and 02 = — 9/2ai to obtain the symmetry 
operator L, ([9]), where ui, U2 are Cartesian coordinates and 



B 



A 




+ ( —ui'^ui u}'^a2Ui H oj^uiu^ f27u? — IGun) | 

V49152 4096 1024 ^ 2;y 2 

/ 59 2 1 2 1 4 4\ 

+ [mE^'' - 4096" "^^^ - 64^" ^^^^J 




C 




-uj'^a2U2 H uj'^U2 : — U2 (243tif — hlQu^u^ + 2)2un) | Lo 

36864 196608 18432 ^ ^ ^ ^ 2; y 2 
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+ UJ Uo UJ OioUo UJ Ur, H 

V 49152 36864 576 V 



+ y-^ulu2U}^a2 + ^uj^uluf - -^U2 (-16^2 + 31ui) uj^j H + (^-^uj^a2ulu2 
- ^uj^ai4 - 1^14 {^4 - '^4) - ^n2 (I14n? - 2<d4) lo""^ L2 

+ -r^u;^U2 U04 + 837uf - 1440^2^?) - 18u^44 - ^^^^2, 
Izo 4U96 



( -^4 - -^4) Ll+u^f -^4 - -^4) HLl - ^_a;2 2^2^2 
V 65536 ^ 49152 V ^ V 16384 ^ 49152 V 12288 ^ 



+ ^^( ^44 - _ i^^A ^2 ^^4^2 (189^2 _ 53^2) 

V2048 2 ^ 144 2 4096 V 1152 ^ 27 ^ 

13 4 4 rT2 ,1a I 113 2 3 2\ 21 4 2 



4 iTT2 if 4 /113 2 2\ 
-UJ Unli + CJ 02 

> 2 V V1024 ^ 256 V 

- (84 - IM44 + 1897xf) + (-A- ^4] u:^\ L, 
32 2 V 2 21 \^64 2 4096 V / 

+ ^4 / iL^2 _ + ia;6nt (I3u2 - 4) + u:^ I ^4 - -^4]] H 

V V256 ^ 256 V 4 2V i 2; \i{)2A ^ 1024 77 

gl3 4 9c 2 / 2 2\ f 4 2025 22 65 4 

- w aiy ^2 - —UJ a2Ui (17ni - IOU2) + uj \ — Ui - ^^%^2 + ^""2 

- 3n?w^u^ (39n? - lOu^) . 



Using Maple, we have checked explicitly that the operator L commutes with H. Note that 
it is of 6th order. Taking the formal adjoint L*, jl9j . we see that S = \{L + L*) is a 6th order, 
formally self-adjoint symmetry operator. 



3 An alternate proof of superintegrability 
for the caged quantum oscillator 

This second proof is very special for the oscillator and exploits the fact that separation of 
variables in Cartesian coordinates is allowed Here we write the Hamiltonian in the form 

Tj o2 I o2 2 2 2 2,4 "1,4 2 

H = d^ + d -ii^x - fi2y + 2 + 2 • 

This is the same as pT|) with ui = x, U2 = y, 4 = —p^uj"^, 4 = —q^uj"^, ai = | — af, 
a2 = \ — 4- look for eigenfunctions for the equation H'i/ = A^* of the form ^' = XY. We 
find the normalized solutions 

Xn = e-^'^i^^'^i+lL^H/Uix'), Ym = e--2^'^y"y-^+-2L-^ifi2y^), 

where the L'^{x) are associated Laguerre polynomials [24] . For the corresponding separation 
constants we obtain 

\x = -2/ii(2n + fli + 1), Xy = -2/X2(2m + 02 + 1). 
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The total energy is — A^; — Xy = E. Taking ni = pfj, and fi2 = qfJ- where p and q are integers we 
find that the total energy is 

E = -2fi{pn + qm+pai + p + qa2 + q)- 

Therefore, in order that E remain fixed we can admit values of integers m and n such that 
pn + qm is a constant. One possibility that suggests itself is that n^n + q, m^m — p. To 
see that this is achievable via differential operators we need only consider 

We now note the recurrence formulas for Laguerre polynomials viz. 



a^^^p (x) = pL^ix) - + a)L°„i(x) = {p + l)L^+i(x) - + 1 + « - x)L^{x) 



Because of the separation equations we can associate p with a differential operator for both m 
and n in the expression for We do not change the coefficients of Lp_^^{x). Therefore we 
can raise or lower the indices m and n in ^ using differential operators. In particular we can 
perform the transformation n ^ n + q, m ^ m — p and preserve the energy eigenvalue. To see 
how this works we observe the formulas 

D+{ni,x)Xn = (^dl - 2xfiA - m + + ^-^^^ Xn = -4/xi(n + l)X„+i, 
D~{i-i2,y)Ym = + 2yfi2dy + /i2 + fJ-W + ^ ^2"^ ^ = -4/i2("i + a2)Ym-l. 

In particular, if we make the choice /xi = 2/i, /i2 = then the operator D~^{2fi,x)D~ {n,y)'^ 
transforms X^Ym to — 128/i^(n+l)(m+a2)(?7i— l+a2)^n+i^m-2- We see that this preserves the 
energy eigenspace. Thus can easily be extended to the case when = p^ and //2 = (ZA* for p and q 
integers. A suitable operator is D~^{piJ,, xYD~ {qn, yY. Hence we have constructed a differential 
operator that commutes with H\ The caged oscillator is quantum superintegrable. This works 
to prove superintegrability in all dimensions as we need only take coordinates pairwise. 

Remark. This second proof of superintegrability, using differential recurrence relations for 
Laguerre polynomials is much more transparent than the canonical operator proof, and it gen- 
eralizes immediately to all dimensions. Unfortunately, this oscillator system is very simple and 
the recurrence relation approach is more difficult to implement for more complicated potentials. 
Special function recurrence relations have to be worked out. Also, we only verified explicitly 
the commutivity on an eigenbasis for a bound state. 



3.1 A proof of superintegrability for a deformed Kepler— Coulomb system 

In [16] there is introduced a new family of Hamiltonians with a deformed Kepler-Coulomb 
potential dependent on an indexing parameter k which is shown to be related to the TTW 
oscillator system system via coupling constant metamorphosis. The authors showed that this 
system is classically superintegrable for all rational k. Here we demonstrate that this system 
is also quantum superintegrable. The proof follows easily from the canonical equations for the 
system. 

The quantum TTW system is H'^ = E'^ with H given by d?]) where 



ui = R, U2 = 9, fi = e ^, f2 = 0, vi = ae 



(21) 
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/3 7 _ 2(7 + /3) 2(7 - /3) cos(2A;g) 

cos2(A:0) sm2(A;0) sin2(2A;^) sin2(2A:0) 

(Setting r = we get the usual expression for this system in polar coordinates r, 9.) In our 
paper [19] we used the canonical form for symmetry operators to establish the superintegrability 
of this system. Our procedure was, based on the results of [20] for the classical case, to postulate 
expansions of F, G in finite series 

F = ^ Aa,b,cEa,b,c{R, 0), G = ^ Bafi,cEa,b,c{R, O), 

a,b,c a,b,c 

Ea^bfi = e^"^ sin''{2ke), Ea,b,i = e^"^ sin^(2fc0) cos{2ke). 

The sum is taken over terms of the form a = + m, b = + n, and c = 0, 1, where m, n are 
integers, oq is a positive integer and bo is a negative integer. Here F, G are the solutions of the 
canonical form equations (fT4|) . (fTSj) . (fT6|) for the TTW system. We succeeded in finding finite 
series solutions for all rational k, and this proved superintegrability. 

In [16] the authors point out that under the Stackel transformation determined by the po- 
tential U = e^^ the TTW system transforms to the equivalent system H"^ = —a^, where 

e^" \ cos^(Kt^) sm^{k6)/ 

Then, setting r = e^^, (/) = 29, we find the deformed Kepler-Coulomb system. 

d -\ 1 di 1 1 ^ = * 

^r^^^^2</' 4r 4r2cos2(A:(/)/2) 4r2sin2(A;0/2)y 4 ' 

From the canonical equations, it is virtually immediate that system (|22p is superintegrable. 
Indeed the only difference between the functions (|2ip defining the TTW system and the functions 
defining the Stackel transformed system is that /i and vi are replaced by fi = exp(4i?) and 
vi = —E exp{2R) (we can set E = H) and the former H is replaced by —a. From this we find 
that the only difference between the canonical equations for the TTW system and the canonical 
equations for the transformed system is that a becomes —H and H becomes —a. Since our 
proof of the superintegrability of the TTW system did not depend on the values of H and a, 
the same argument shows the superintegrability of the transformed system. Also, any solution 
(i.e., a second commuting operator) for the TTW system gives rise to a corresponding operator 
for the transformed system by replacing a and H with —H and —a, respectively. 

Example. For the deformed Kepler system with k = 2, we take p = 2 and q = 1 and use 
Gaussian elimination to find a solution with nonzero coefficients 

^-2,1,0, ^-1,1,0, -B_2,o,o> -B_2,o,i, -B-1,0,0, -B-1,0,1, -Bo,o,i, 

in which A_2,i,o and A-i^i^ are two independent parameters. To achieve the lowest order 
symmetry and ensure that the A, B, G and D are polynomial in L2 and H, we choose A-2,1,0 = 
32(^2 — 4) and A^i^i^ = and find, 

A = 32e-^^ sin(46')L2 - 128e"''^ sin(46l), 

B = 8e-^^ cos{49)Ll + (8(/3 - 7)6"^^ + AEe'^^ cos(40) - Se"^^ cos(40))L2 
- 96e~^^ cos(40) - leEe'^^ cos(40) + 40(7 - I3)e-^^ + 4E{/3 - -/)e-^^, 

C = -8e-^^sm{49)Ll - 4sin(46l)i?L2 + 8(6"^^ - Ee''^^) sm{49)L2 + 20sin(46l)i? 
+ 96e-^^ sin(46l) + 32Ee-^^ sm{49) - E^ sin(46l). 
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D = -32e~'^^cos(46l)L^ -8cos{49)HL2 

+ (128e-^-^ cos(40) - UEe'^^ cos{Ae) + 16(7 - /3)e-^^)L2 + 2Acos{Ae)H 
+ 48(/3 - 7)e~^^ + A8Ee-^^ cos(40) - 4£;(7 - P)e-^^ + 2E^ cos(40). 

It has been explicitly checked, using Maple, that the fifth order operator obtained from these 
expressions commutes with H. 

3.2 A proof of superintegrability for the caged oscillator on the hyperboloid 

Using this same idea we can find a new superintegrable system on the 2-sheet hyperboloid by 
taking a Stackel transformation of the quantum caged oscillator = E^, ()17p in Cartesian 
coordinates ui, U2 by multiplying with the potential U = ^/u\. The transformed system is 



ul {dl + dl + uo'^ulip^ul + q^ul) - Eul + 03^^ ^ = -oi^. 



We embed this system as the upper sheet sq > of the 2-sheet hyperboloid Sq — — s| = 1 in 3- 
dimensional Minkowski space with Minkowski metric ds^ = —ds^ + ds^ + ds^, via the coordinate 
transformation 

l + u\+ul 1 - n? - It, U2 

So = , si = , S2 = — , ui> 0. 

lui zui Ul 

Then the potential for the transformed system is 

v=^- + - + 73^^ ■ (23) 

si (S0 + Sl)2 (S0 + Sl)4 (S0 + Si)3 

This is an extension of the complex sphere system [S2], distinct from the system ([6]) that we 
proved classically superintegrable. It is an easy consequence of the results of [T7] that the 
classical version of system (j23p is also superintegrable for all rationally related p, q. However, 
quantum superintegrability isn't obvious. However, from the results of Section [2.21 it is virtually 
immediate that this new system is quantum superintegrable for all relatively prime positive 
integers p, q. This follows from writing down the canonical equations (jlSp . (jl6p . first for the 
caged oscillator where 

-PI f n 222, "1 222, 0^2 

and then for the Stackel transformed system where 

f ^ f n 222771 222, '^2 

71 = —, 72 = 0, Vi=ujpu;^-E, V2=ujqu2 + —. 



u 



The equations are identical except for the switches ai — )■ —E, H — )• —ai. Thus our proof of 
quantum superintegrability for the caged oscillator carries over to show that the system on the 
hyperboloid is also quantum superintegrable. 



4 Discussion 



A basic issue in discovering and verifying higher order superintegrabily of classical and quantum 
systems is the difficulty of manipulating high order constants of the motion and, particularly, 
higher order partial differential operators. We have described several approaches to simplify such 
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calculations. Although our primary emphasis in this paper was to develop tools for verifying 
classical and quantum superintegrabity at all orders, we have presented many new results. The 
classical Eucidean systems [E8], [E17] and most of the examples of superintegrability for spaces 
with non-zero scalar curvature are new. We have explored the limits of the construction of 
classical superintegrable systems via the methods of Section 11.11 for the case n = 2. Further 
use of this method will require looking at n > 2, where new types of behavior occur, such as 
appearance of superintegrable systems that are not conformally flat. We also developed the 
use of the canonical form for a symmetry operator to prove quantum superintegrability. We 
applied the method to the n = 2 caged anisotropic oscillator to give the first proof of quantum 
superintegrability for all rational k. We used the Stackel transform together with the canonical 
equations to give the first proofs of quantum superintegrability for all rational A: of a 2D deformed 
Kepler-Coulomb system, and of the caged anisotropic oscillator on the 2-hyperboloid. Then 
we introduced a new approach to proving quantum superintegrability via recurrence relations 
obeyed by the energy eigenfunctions of a quantum system, and gave an alternate proof of the 
superintegrabilty for all rational k of the caged anisotropic oscillator. This proof clearly extends 
to all n. The recurrences for the caged oscillator are particularly simple, but the method we 
presented shows great promise for broader application. Clearly, we are just at the beginning 
of the process of discovery and classification of higher order superintegrable systems in all 
dimensions n > 2. 
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